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£> ' Abstract 

PsJ ■ We consider the production of strange particles in Pb-Pb and p-A collisions at 

the SPS energy reported by the WA97 experiment. We show that the observed 
enhancement of strange baryon and antibaryon yields in Pb-Pb collisions relative 
to p-Be and p-Pb can be explained in terms of the statistical model formulated in 
canonical ensemble with respect to strangeness conservation. The importance and 
' the role of strangeness under saturation is also discussed. 

a: 

1 Introduction 

Various experiments at CERN SPS reported an enhancement of strange particle produc- 
tion when going from p-p, p-A to A-A collisions [1], |2], In particular the recent 
analysis of Pb-Pb reaction by the WA97 collaboration has shown, by comparing p-Be, 
p-Pb with Pb-Pb, a strong enhancement in the pertinent strange and antistrange baryon 
yields [0, 0, §]• In addition this enhancement was found to increase with strangeness 
content of the particle. 

In this letter we show that the observed features of the data can be well described 
by the thermal model for particle production using canonical formulation of strangeness 
conservation. 

The statistical models formulated in the grand canonical (GC) ensemble, with respect 
to conservation laws, were shown to give a satisfactory description of particle production in 
heavy ion collisions at SPS and AGS energies || [7|, |9|, [TOfl . However, when applying this 



model to the low energy heavy ion [|TT], or high energy hadron-hadron or e e~ reactions 
one necessarily needs to treat the conservation laws exactly [|TJ| . The exact conservation 
of quantum numbers, that is the canonical (C) approach, is known to reduce the phase 
space available for particle production due to additional constraints appearing through 



requirements of local quantum number conservation |T3|, [TJ|, |J| [16! ■ Thus, the immediate 
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question is to what extend the canonical suppression can account for the observed effect 
in the WA97 experiment. 

The concept of canonical suppression as an origin of strangeness enhancement has 
been already addressed in the literature W7\ and also tested in the context of the SPS 



data UlSfl . Indeed, it was shown |18| that the increase of Kg/h~ and (A + K)/tt ratios 
from p-p to A- A collisions can be possibly explained by the canonical effect. In view of 
the new data from the WA97 collaboration we will test further this idea and show that 
the suppression of strange particles phase space when going from GC to C formulation of 
strangeness conservation can quantitatively explain the reported data. 

In the following section we present the main aspects of the statistical model formulated 
in canonical ensemble and then we compare the model with the experimental data. 



2 Canonical Formulation of Strangeness Conserva- 
tion 

We consider a gas composed of all particles and resonances with mass up to m ~ 2.4GeV. 
Conservation of baryon number is introduced in GC ensemble, thus it is controlled by 
baryochemical potential. Strangeness conservation, however, is introduced on the canon- 
ical level. Having in mind the application of the model to particle production in heavy 
ion collisions we also require that the overall strangeness of the system is zero. Using 



the general procedure for the canonical strangeness conservation [15, 16] the partition 



function of a gas of particles with strangeness 0, ±1, ±2, ±3 can be written as follows: 



Zo=^~ [* d<j> exp(]T S s e is % (1) 

27T J 7r s= _ 3 



where S s = Vq J2i an d the sum is taken over all particles and resonances carrying 
strangeness s. For particle of mass rrii, degeneracy factor dj, baryon number Bi, and 
assuming Boltzmann statistics we have: 

dl -m\TK 2 (rrii/T) exp (B if i/T) . (2) 
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With the particular form of the partition function given by Eq. ([]]) the density n s of 
particle i with strangeness s in volume V is obtained by the replacement i— > AjZ, in 
Eq. (|IJ) and then taking an appropriate derivative [T(| : 



(N ±s ) \j dhiZp , (S ±1 

n ±s = t r = TT^T, A 4 =l - Z ±s 



where 



v v o\ IAi " 1 ^ s (s +1 s^y/ 2 
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x k = 2\/S k S„ k , A = 7 ^, (4) 
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Figure 1: Volume dependence of strange baryon densities normalized to their value 
calculated in V = V = 7.4 fm 3 for T = 168 MeV and fi = 266 MeV. 

and the partition function 

oo 

Z ~ I (xi)I (x 2 ) + X) I 2m (a:i)/ m (a; 2 )[A m / 2 + A~ m,2 \. (5) 

m=l 

In the derivation of Eq. (|3]) we have neglected, after differentiation over particle fu- 
gacity, the term S±^. This approximation, however, due to small value of S±z coefficients 
in comparison with S±\ and S± 2 is quite satisfactory. The complete expression for the 



partition function without any approximation can be found in reference [|L6] . 

Due to the locality of strangeness conservation, we have introduced two volume pa- 
rameters, V and Vo in Eq. (|3|) fL3f . The parameter V describes the volume of the system 
at freeze out, thus V is also the overall normalization factor which fixes the particle mul- 
tiplicity from the corresponding density. The volume Vo is interpreted as a correlation 
volume where particles are created to fulfill the requirement of local conservation of quan- 



tum numbers. It was shown Jy]] that in the canonical description of particle production 
in low energy heavy ion collisions V ^ Vq is required to reproduce both the particle mul- 
tiplicity ratios and the particle yields. The volume parameter V scales with the total 
number of participants, V ~ A part , whereas V = 47i\Rq/3 where R ~ (1.15 ± 0.05)A^ 3 
with A p being the number of the projectile participants. 

Comparing the result of Eq. @ with the standard expression for the particle density 
in GC ensemble, one can see that the term multiplying Z± s in Eq. (H) plays the role of 



the strangeness fugacity |TT], [L5j. Indeed taking the limit for large x iy in Eq. (^|), that is 



the limit of high temperature and/or large Vq, the GC result for the particle density is 
recovered. In the opposite limit the C ensemble leads to a strong suppression of strange 
particle densities. 

In Fig.|l] we show the dependence of the density of strange particles on Vq for fixed 
T = 168 MeV and /i = 266 MeV. In order to illustrate the enhancement with increasing 
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volume and strangeness content of the particle we show in Fig.[l] the density of A, S~ and Q 
baryons normalized to its value calculated in the volume Vq = Vq = 7.4 fm 3 , as a function 
of Vo/Vq. It is clear from Fig.|l| that the particle densities are increasing with the size of the 
system. For sufficiently large Vo the densities are saturating which indicates that the GC 
limit is reached. The enhancement is seen to increase with the strangeness content of the 
particle and obviously depends on the choice of Vq. Keeping the value of baryon chemical 
potential and/or temperature as being Vq independent, the enhancement of antibaryons is 
similar to the one seen in Fig.|l| for baryons. Thus, the possible asymmetry in enhancement 
for strange baryons and antibaryons when going from small to large system requires fi 
and/or T to be Vo dependent. 

In the following section we will discuss to what extend the above model can be used 
to understand strangeness enhancement versus centrality for baryon and antibaryons re- 
ported by the WA97 collaboration. 

3 Comparison with the Experimental Data 

It is seen in Fig.|l] that the model reproduce the main features of centrality dependence of 
strangeness enhancement in the WA97 data. The quantitative comparison of the model 
with the experimental data require, however, a particular choice of the parameters, that 
is: the temperature T, the volume Vo and the baryon chemical potential /i. We have used 
for Pb-Pb collisions: T = 168 MeV and fi = 266 MeV, the values which have been shown 
|| to give a good description of all data on particle multiplicity ratios measured by the 
NA49 and WA97 experiments. In general temperature and chemical potential could vary 
with Vo that is also with A p . From previous analysis we know ||19| that for the same 
collision energy the temperature variations with A part is small. Thus, we keep T to be the 
same for both, p-Be and Pb-Pb systems. The chemical potential, however, as discussed 
in the previous section, has to be taken as A part dependent since the data are showing 
larger enhancement for strange than for antistrange baryons. This indicates lower value 
of /j, in p-Be than in Pb-Pb system. We use Jx in p-Be as a free parameter which is chosen 
to reproduce the enhancement asymmetry between S~ and S + or simply the ratio 
measured in p-Be collisions. 

Following the interpretation of V , described in the previous section, we chose R ~ 1.2 
fm in p-Be that is Vq = Vq. To make a quantitative comparison of the model with the 
experimental data we are left with only one parameter, the effective chemical potential Ji 
in p-Be which was taken to be 150 MeV. This value is required in comparison with /i = 266 
MeV in Pb-Pb in order to obtain the correct splitting between baryons and antibaryons 
of the enhancement factors without affecting the average enhancement factors. 

In Fig.|2| we show the densities calculated in C ensemble for large Vo normalized to 
its value at Vo in to p-Be system. This corresponds to the experimental comparison of 
yield/participant of Pb-Pb to p-Be system. The full lines in FigPJ are calculated from 
Eq. (|3p by varying Vq and keeping /i = 266 MeV. The change of Vq with centrality in Pb- 
Pb collisions is included relating Vq = 47ri?Q/3 with A part through, R ~ 1.2(y4 part /2) 1 / 3 , 
which determines approximately the number of projectile participants. 

The actual path from grand canonical to canonical result for the particle densities 
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Figure 2: Particle yields/participant for strange baryons and antibaryons calculated in 
the canonical ensemble with Vq ~ A part /2 normalized to their values at Vq. Also shown 
are the WA97 experimental data for yields per participant in Pb-Pb relative to p-Be 0. 
For description of curves see text. 



would required an additional knowledge of A part dependence of chemical potential which is 
effectively changing from 266 MeV in Pb-Pb to 150 MeV in p-Be system. Formulating, in 
addition to strangeness, also baryon number conservation canonically, the path from GC to 
C result as a function of Vq can be established [20|. In the simplified model, presented here, 



where baryon number conservation was introduced effectively trough chemical potential 
we can only study the relative enhancement of strange baryon densities from small p-Be 
to large Pb-Pb system. 

The result presented in Fig.|2] shows that the relative enhancement of strange baryon 
and antibaryon yield/participant from p-Be to Pb-Pb, found experimentally, can be well 
described by the effect of canonical suppression and the change of \x from 266 MeV in 
Pb-Pb to 150 MeV in p-Be collisions. The enhancement appears in the model to be 
almost saturated for A part ~ 100 which indicates that the system is close to its GC limit. 
The enhancement pattern for different strange baryons and antibaryons densities is also 
understood within the model as a consequence of different dependence on Vq. In the large 
system like Pb-Pb, and for large collision energy required to maintain high T, the density 
n s of particle carrying strangeness s, is Vq independent. In the opposite limit, however, 
this dependence is changed to n s ~ V S which can be verified from Eq. ([|). For small Xi 
we have approximately: 



J ±s 



(S 



Expanding the Bessel functions I n (x) 



(s +1 s^y/ 2 i ( Xl y 

x n in Eq. @ we see that n s 



(6) 
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The variation of Vq with the number of projectile participants and the assumption that 
V 7^ Vq is also justified by the data. Indeed the recent data of the WA97 collaboration 
show that strange particle yield/v4 part is independent on A part in p-A collisions. This is 
well reproduced by canonical treatment since Vq is the same for all A in p-A collisions. 
If one could, however, take V = Vq then, since V ~ A part , a strong variation with A part 
of strange particle densities should be observed in p-A collisions when changing A. In 
addition this variation would depends on the strangeness content of the particle which 
is not experimentally observed in p-A. Thus from the point of view of the considered 
canonical model, p-Be and p-Pb are the same system with the only difference being 
contained in slightly larger Jx in p-Pb then p-Be. We found that the chemical potential in 
p-Pb increases by about 10% with respect to p-Be. This explains similar excess of particle 
yield/ A par t in Pb-Pb relative to p-Be and p-Pb collisions. 

Until now we have assumed that p-A as well as Pb-Pb system can be described by 
the statistical model in chemical equilibrium. For Pb-Pb collisions chemical equilibrium 
assumption has been found to give a good agreement of the model with the experimental 
data 0. In p-A collisions, however, due to small time scale, deviations from strangeness 
chemical saturation cannot be excluded. In a previous analysis |12| of particle production 
in p — p, p — p and e + — e~ collisions in terms of C ensemble strangeness under saturation 
was accounted for by introducing an additional parameter 7 which measures deviation of 
strangeness from full equilibrium pl |. Also in this approach Vq was assumed to coincide 



with V. In the following we discuss how the possible out of chemical equilibrium together 
with the V — Vq ~ A part assumption will modify our results and conclusions. 

To introduce non-equilibrium effects via 7 parameter is quite straightforward |ZT], [l(| . 
It amounts to the replacement S± s 1— > , y s S± s in Eq. (|3p. First we observe that applying the 
limit X\ 1 — ^ in Eq. (||) the density n s of particle with strangeness s, scales as {^Vy . For 
larger Xi we found that, in the parameter range relevant for p-A system, this dependence 
is changed approximately to n s ~ {^ 2 - 7 V) s . Thus, keeping 7 2 - 7 V r = Vo, the value used in 
equilibrium model, we can as well reproduce the experimental data shown in Fig.^| in non- 
equilibrium formulation. In this approach, however, in order to reproduce experimentally 
observed linear dependence of strange particle yields with A part in p-A collisions one 
needs, for each value of A part , to adjust 7 such that 7 2,7 A part = const. Obviously, this 
requires a decreasing 7 with A part . In general, one would rather expect the opposite 
behavior. In order to ovoid, however, the possible ambiguity with respect to strangeness 
undersaturation factor 7 one needs to perform an additional detailed analysis in proton 
induced collisions for different targets and collision energies. This analysis could possibly 
settle to what extend there is a need for genuine global strangeness enhancement related 
with the change in 7 between p-A and A-B collisions. 



4 Summary 

We have used the comprehensive set of the recent results of the WA97 collaboration on 
the production of strange and multistrange hadrons in p-A and Pb-Pb collisions in order 
to check the validity of statistical model description of strangeness enhancement central- 
ity dependence. The statistical model was verified to give a good agreement with all 
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experimental data in A-A collisions at SPS energies. From the point of view of statistical 
mechanics the application of the model for smaller system like p-A requires obvious modi- 
fications. In the large system strangeness is conserved on the average, thus it is controlled 
by strange chemical potential in the grand canonical ensemble. In the small system, 
however, strangeness conservation laws must strictly be obeyed in each single collision 
which requires canonical ensemble. We have shown that the required passage from grand 
canonical to canonical ensemble reproduces the basic features of strangeness suppression 
from Pb-Pb to p-Be or p-Pb collisions reported by WA97 collaboration. The observed 
asymmetry in the enhancement of strange baryons and antibaryons requires additional 
assumption on A part dependence of baryon chemical potential /i. We have shown that 
with the particular choice of /i asymmetry in p-A and Pb-Pb systems the model provides 
good quantitative description of the data. 

Deviation from chemical equilibrium has been included in our discussion only in p- 
A system. Recently an interesting non-equilibrium model for the description of particle 
production in A-B collisions has been proposed f22|| . The interpretation of centrality 
dependence of strangeness enhancement in a consistent non-equilibrium picture would 
require additional analysis not presented in this paper. 
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